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(U)  It  is  well  known  that  the  static  pressure  coefficient  values  obtained 
from  tests  of  axisymmetric  bodies  in,  for  example,  a large  water  tunnel 
indicate  the  existence  of  tunnel  wall  interference.  This  interference  is 
due  to  blockage  experienced  by  the  body  operating  within  the  boundaries  of 
the  test  section  walls,  and  to  the  skin  friction  on  both  the  tunnel  walls 
and  the  surface  of  the  body,  which  in  turn  causes  finite  thickness  boundary 
layers  to  develop.  These  factors  result  in  a decrease  in  static  pressure  . : 
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along  the  test  section  which  leads  to  a spurious  horizontal  buoyancy  on  the 
body. -^These  interferences  become  more  important_  as  J:he  size  of  the  body 
increases  relative  to  the  fixed  dimensions" ~dF~*the  test  sections*  Potential 
flow  calculations  can  only  account  for  solid  blockage  effects.  A control 
volume  analysis  has  thus  been  carried  out  from  which  pressure  coefficient 
correction  terms  due  to  horizontal  buoyancy  have  been  obtained.  They  may 
be  directly  applied  to  any  inviscid  pressure  distribution  calculated  for 
any  given  body  operating  in  a given  tunnel.  These  correction  factors  are 
computed  for  a large,  streamlined  body  operating  over  a range  of  velocities, -i 
in  The  Pennsylvania  State  University,  Garfield  Thomas  1.22-m.  diameter 
Water  Tunnel.^ They  are  applied  to  a theoretically-determined  potential  flow, 
in-tunnel  body  pressure  distribution  to  obtain  the  viscous  flow  pressure 

distributions  which  are  compared  to  those  determined  experimentally.  ^The  ~ 

agreement  between  theory  and  experiment  is  shown  to  be  very  good  over  all 
speeds  considered.  Traditional  methods  of  accounting  for  horizontal  buoyancy 
which  are  obtained  from  measurements  of  the  longitudinal  static  pressure 
gradient  in  a clear  tunnel,  are  applicable  only  to  small  bodies.  The  present 
method  permits  one  to  treat  relatively  large  bodies  which  have  a much  stronger 
interaction  with  the  tunnel  walls. 
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Abstract:  It  is  well  known  that  the  static  pressure  coefficient  values 
obtained  from  tests  of  axisymmetric  bodies  in,  for  example,  a 
large  water  tunnel  indicate  the  existence  of  tunnel  wall  inter- 
ference. This  interference  is  due  to  blockage  experienced  by 
the  body  operating  within  the  boundaries  of  the  test  section 
walls,  and  to  the  skin  friction  on  both  the  tunnel  walls  and 
the  surface  of  the  body,  which  in  turn  causes  finite  thickness 
boundary  layers  to  develop.  These  factors  result  in  a decrease 
in  static  pressure  along  the  test  section  which  leads  to  a 
spurious  horizontal  buoyancy  on  the  body.  These  interferences 
become  more  important  as  the  size  of  the  body  increases  relative 
to  the  fixed  dimensions  of  the  test  section.  Potential  flow 
calculations  can  only  account  for  solid  blockage  effects.  A 
control  volume  analysis  has  thus  been  carried  out  from  which 
pressure  coefficient  correction  terms  due  to  horizontal  buoyancy 
have  beer  obtained.  They  may  be  directly  applied  to  any  inviscid 
pressure  distribution  calculated  for  any  given  body  operating  in  a 
given  tunnel.  These  correction  factors  are  computed  for  a large, 
streamlined  body  operating  over  a range  of  velocities  in  The 
Pennsylvania  State  University,  Garfield  Thomas  1.22-m.  diameter 
Water  Tunnel.  They  are  applied  to  a theoretically-determined 
potential  flow,  in-tunnel  body  pressure  distribution  to  obtain 
the  viscous  flow  pressure  distributions  which  are  compared  to 
those  determined  experimentally.  The  agreement  between  theory 
and  experiment  is  shown  to  be  very  good  over  all  speeds  considered. 
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Nomenclature 

A the  blocked  area  between  the  body  and  duct  wall 

Ag  the  effective  area  between  the  body  and  duct  wall  including 
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A area  of  duct 

o 
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£ distance  from  body  nose  to  D 

n J max 

Pq  free-stream  static  pressure 
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Pb 

K 

\ 


u 

oo 


local  fluid  pressure 

pressure  exerted  by  the  body  on  the  fluid 

"displaced-body"  pressure  exerted  by  the  body  on  the  fluid 

radius  of  curvature  of  the  body  tailcone 

body-length  Reynolds  number 

radius  coordinate 

arc  length  coordinate 

transition  arc  length 

free-stream  velocity  in  duct  far  upstream  of  body 
local  fluid  velocity 

velocity  profile  for  body  (i  = 1,  laminar;  i = 2,  turbulent) 

reference  velocity  in  potential  flow  computation 

numerical  rake  output  far  upstream  of  body  (U^/u^) 

numerical  rake  output  near  the  body  at,  for  example,  a pitot 
probe  location 

axial  coordinate  where  origin  is  i upstream  from  nose  of  the  body 

body  axial  coordinate  (x  - Jt^) 

non-dimensional  body  axial  coordinate  (x/&n) 

dummy  variable  of  integration 

non-dimensional  body  radius  (r^/Dmay/2) 

boundary-layer  thickness 

body  boundary-layer  thickness  (i  = 1,  laminar;  i = 2,  turbulent) 
boundary-layer  displacement  thickness 
the  expression  given  by  Eq.  (47) 
kinematic  viscosity 
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p fluid  mass  density 

T wall  shear  stress 

<J>  angle  formed  between  tangent  line  of  body  and  the  x-axis 

Special  Subscripts 


b 

refers 

to 

body 

E 

refers 

to 

experimental 

I 

refers 

to 

inviscid 

m 

refers 

to 

measured 

V 

refers 

to 

viscous 

w 

refers 

to 

duct  wall 

Ml. 
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INTRODUCTION 

The  investigation  of  flows  about  bodies  of  revolution  in  cylindrical  ducts 
is  of  general  interest  to  those  concerned  with  turbomachinery,  wind  tunnel,  and 
water  tunnel  testing,  to  note  a few.  There  have  been  several  methods  proposed 
for  the  determination  of  the  potential  flow  about  bodies  operated  within  ducts 
of  both  consta.it  and  non-uniform  cross  section  [1-5].  These  methods  allow  for 
blockage  effects  due  to  the  duct  (or  test  section)  walls  to  be  accounted  for. 
However,  these  methods  do  not  allow  for  viscous  effects  that  occur  on  the  body 
and  tunnel  walls.  Viscosity  gives  rise  to  skin  friction  on  all  surfaces  in  con- 
tact with  the  flowing  fluid;  which,  in  turn  gives  rise  to  boundary  layer  growth. 

The  displacement  thickness  varies  from  the  beginning  of  boundary 
layer  growth  to  distant  points  far  downstream, effectively  changing  the  cross- 
sectional  areas  of  both  the  body  and  the  duct.  This  added  constriction  of  the 
flow  and  the  skin  friction  creates  a second  type  of  blockage  characterized  by  a 
static  pressure  drop  from  the  beginning  of  boundary  layer  growth  to  any  arbitrary 
point  downstream.  This  net  change  in  static  pressure  may  be  interpreted  in  the 
form  of  a force  in  the  horizontal  direction  (assuming  the  tunnel  test  section 
is  horizontal)  which  is  the  so-called  "horizontal  buoyancy." 

In  many  applications,  it  is  necessary  to  know  as  precisely  as  possible,  the 
pressure  distribution  on  a particular  body  operating  in  a given  tunnel  test  sec- 
tion. Because  of  horizontal  buoyancy  effects,  a typical  inviscid  flow  calculation 
will  not  yield  totally  accurate  results  for  relatively  large  diameter  bodies  (rela- 
tive to  the  diameter  of  the  test  section) . This  statement  may  be  verified  through 
a direct  comparison  of  experimentally  determined  pressure  distributions  with  a 
potential  flow  prediction.  Figure  1 shows  such  a comparison.  The  axisymmetric 
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body  here  is  3.05-m.  long  and  has  a maximum  diameter  of  32.39-cm.  The  . • 
pressure  coefficient  data  were  obtained  on  this  body  while  it  operated  in 

V V 

the  Garfield  Thomas  48-inch  (1.22-m.)  diameter  Water  Tunnel  lcated  at  the  Applied 
Research  Laboratory  of  The  Pennsylvania  State  University.  A description  of 
the  experiments  will  be  presented  in  a later  section.  The  potential  flow 
calculation  shown  in  Figure  1 was  performed  using  the  Douglas-Neumann  code 
[1,2^.  The  solid  blockage  due  to  the  1.22-m.  diameter  test  section  is  included 
in  this  computation.  The  experimental  pressure  coefficients  are  seen  to  be 
lower  than  those  predicted  at  a given  point  on  the  body  and  are  velocity  - 
dependent.  This  discrepancy  is  due  to  horizontal  buoyancy. 

Measurements  of  body  pressure  distributions  are  not  economically  practical 

'1 

in  all  cases.  The  next  best  alternative,  therefore  is  to  approach  analytically 
the  problem  of  horizontal  buoyancy.  This  report  presents  such  an  analysis,  where 
. correction  terms  are  derived  that  may  be  applied  to  the  results  of  a given 

potential  flow  calculation.  As  examples  of  the  analysis,  we  include  calcula- 
tions for  the  b ody  of  Figure  1 and  compare  them  with  the  experimentally 
determined  data.  Although  the  analysis  given  here  is  for  axisymmetric  bodies 
operating  in  cylindrical  test  sections  of  constant  radius,  it  is  quite  general 
so  that  its  adaptation  to  other  configurations  (variable  radius  test  sections, 
for  example)  would  seem  straightforward. 

ANALYSIS 

Potential  Flew  Considerations 


The  potential  flew  pressure  distribution  on  an  arbitrarily-  shaped  axisymmetric  body 
operating  within  the  confines  of  a cylindrical  duct  of  circular  cross  section  may  be 
calculated  using  the  Douglas-  Neumann  computer  code  [1,2],  among  others  [3-5]. 
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A typical  computation  with  the  Douglas-Neumann  code  might  use  the  modeling  shown 

in  Figure  2(a).  Because  of  the  limited  core  size  on  our  digital  computers,  the 

length  of  the  duct  (which  is  described  by  discrete  points)  must  be  finite  when 

using  this  method.  If  the  duct  is  assumed  to  extend  from  x=-£.tox*=2  + L, 

where  L is  the  body  length,  then  the  resulting  pressure  coefficient  values 

on  the  body  (Cp)  are  normalized  to  an  unspecified  reference  velocity,  u . We 

o 

know,  however,  that  in  the  limit  2 -*■  <»,  uq  -*•  (UB  is  the  true  velocity  at 
infinity).  Consequently,  in  conjunction  with  our  computation  of  the  potential 
flow  about  the  body,  we  also  determine  ur  = U at  x * -t,  where  2 is  assumed 
much  larger  than  L.  Knowing  u^,  the  calculated  pressure  coefficient  values  may 
then  be  corrected  to  the  reference  velocity  by  a simple  computation: 


Cp  = 1 - 


' 

UT 

>2 

• * 

u 

I 

o 

u 

u 

o 

00 

> 4 

1 4 

1 - (1  - Cp) / u 


(1) 


where  Cp^.  represents  the  potential  flow  pressure  coefficient  on  a body  operating 
in  a cylindrical  tube  of  infinite  extent.  Equation  (1)  is  most  useful  to  use 
when  a number  of  pressure  distributions  are  to  be  calculated  for  a family  of 
bodies  and  compared;  it  assures  us  of  a consistent  reference  velocity. 

By  defining  a control  volume  from  Figure  2(a)  to  extend  from  x « -®  to  an 
arbitrary  x on  the  body,  the  continuity  equation  may  be  written  as 


A U 
o 00 


A u„ 


where 


ir  r 


w 


(2) 


(3) 


A “ * (*w1 2  - rb2) 


(4) 
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with  and  being  the  tunnel  wall  and  body  radii,  respectively*  Vfe  may 

also  write  the  momentum  flux  equation: 


p Uj2  A - p U^2  A0  " F!  + f2  + f3  . (5) 

where  p is  the  density  of  the  fluid,  and  the  F^s  (i  = 1,  2,  3)  are  the 
forces  acting  on  the  control  volume.  Assuming  positive  forces  to  act  in 
the  positive  x-direction,  we  have: 


and 


F - P A 
1 o o 


(6) 


F2  “ -Px  A • (7) 

where  Pq  is  the  free-stream  static  pressure,  and  p^  is  the  potential  pressure 
that  occurs  in  the  fluid  at  any  arbitrary  location  within  the  control  volume.  The 
third  force,  is  due  to  the  pressure  exerted  by  the  body  in  the  x-direction. 

It  depends  therefore  on  the  pressure  distribution  of  the  body,  p^(x).  The 
differential  form  of  this  force  may  be  written  as 

dF^  = - p^  sin<j>  dA^  , (8) 

where  <J>  is  the  angle  between  the  tangent  line  of  the  body  at  point  x and  the 
x-axis , and 

= 2tt  r^  ds  = 2tt  r^  dx/cos<J>  , (9) 

where  s is  the  arc-length  coordinate.  We  can  integrate  dF^  along  the  body 
and  get 


F^  = -2it  | p^  r^  tan<f>  dx'  , (10) 

where  x'  is  a dummy  variable  of  integration. 

We  may  now  solve  Eq.  (5)  for  p^  (contained  in  F^)  and  use  Eq.  (2)  to  obtain: 

A A 2 a ^ 

pi " p0  ir  - p u»2  -j- + p u«>2  r - F J pb  rb tan ♦ dx'  • (11) 

o 


IH 


ifinliillii 


where  u is  the  non-uniform  velocity  profile  between  the  body  and  tunnel  wall. 
In  like  manner,  we  can  write  the  change  in  momentum  flux  as 


r 
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where 


P A 
o o 


- Pv  A . 


— 2tt  p'  r,  tan<{>  dx' 


-2"  rb  Tb  dx' 


F8  - -2*  r, 


X 

i T dx' 
w J w 


The  quantity  pv  is  the  viscous  pressure  that  occurs  in  the  fluid  between  the 

body  and  tunnel  wall  at  location  x , and  p£  is  the  "displaced-body"  pressure, 

i.e.,  it  is  a perturbed  due  to  the  added  acceleration  of  the  flow  caused 

by  the  body  boundary-layer  development.  It  may  be  calculated  using  potential 

* * 

flow  methods,  but  with  the  body  radius  re-defined  by  r'  = r,  + 6/,  where  6. 

b b b b 

is  the  body  boundary-layer  displacement  thickness.  In  Eq.'s  (17)  and 
(18) , and  x^  are  the  wall  shear  stresses  that  occur  on  the  body  and  tunnel 

wall,  respectively.  It  is  noted  that  does  not  contain  a cos<J>  term  because 

Fy,  by  assumption,  is  the  x-component  of  the  force  due  to  body  wall  shear  stress. 
Equation  (13)  can  be  readily  solved  for  p^  to  obtain  an  analogous  form  of 
Eq.  (11)  for  viscous  flow. 


Pressure  Coefficient  Correction 


To  obtain  a viscous  flow  correction  for  the  in-tunnel,  inviscid  pressure 
coefficient,  Cp^  , we  form  the  difference 


j 
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x 

pfa  rb  tan<f>  dx'  - j rb  tan<|)  dx' 


x x 

■ J rb  Tb  dx'  - rw  | T„  dx'  - P | 


dx'  - p u rdr 

w 

o r. 


(19) 


where  the  change  of  variable  x = x + l was  made  in  Eq.  (11), 


A pressure  coefficient  correction  term  may  be  defined  by: 


ACp  «=  Cpv  - Cpj 


(20) 


where 


Cp 


PVPo 


V 1 


2 p Uc 


(21) 


and 


Cp 


PI  ~ Po 


I 1 


2 p U» 


(22) 


such  that 


AC 


Pv  - Pi 


P I 


2 PU» 


(23) 


It  is  convenient  to  express  ACp  as  the  sum  of  individual  terms,  i.e., 


ACr 


l ACp 
i-1  1 


(24) 
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where 


ACp.  - 2 
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(27) 


and 


4ir 

ACp, J 

A p 


X X 

p^  tan<{>  dx'  - j p£  tan<J>  dx' 


(28) 


Under  the  assumption  that  a turbulent  boundary  layer  originates  at  x = 0,  and 

that  (SL  + L)  ~ 10  r , then  the  wall  shear  stress,  T , can  be  calculated  from 
o w w 

flat-plate  expressions  [6,  7].  After  making  the  substitution  [6] 


T = 0.0225  p U 
w o 


6 

00  w 


1/4 


(29) 


into  Eq.  (27),  where  the  boundary-layer  thickness  is  given  by 


6 ■ 0.37  x 
w 


U_  x 


V-l/5 


(30) 


we  can  evaluate  the  integral  of  Eq.  (27)  in  closed  form.  Thus,  we  find  that 


0.144  tt 
“*3  ’ ‘r-  — — 


. «L  ) 


1/5 
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(31) 
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where  V is  the  kinematic  viscosity,  and  A * Aq  for  x SLq.  For  arbitrarily- 
shaped  bodies,  the  integrals  of  Eq.'s  (26)  and  (28)  must  be  evaluated  numerically. 

The  integral  of  Eq.  (25),  however,  may  be  evaluated  in  closed  form  if  we  make 
certain  assumptions  regarding  the  velocity  profile  u(r)  . Clearly, 


V*i 


r -6 
w w 


f 

2 f 

2 

2 

f 

5J 

u rdr  = 

Ui 

rdr  + uy 

rdr  + 

-v  r,+6. 

b b b i w w 

where  subscript  i equals  1 for  laminar-body  flow  and  is  equal  to  2 for 

turbulent-body  flow.  The  boundary-layer  velocity  profile,  u , at  the  tunnel 

w 

wall  is  (as  above)  assumed  turbulent,  i.e.,  [6], 

ft  -r]l» 
w 


l s“« 

w w 


rdr  ,(32) 


Uw  = “v  [ ~ 

On  the  turbulent-flow  body,  we  have: 


•\l/7 

r - rb  ' 

u2  = “v 


where  6^  is  of  the  same  analytical  form  as  Eq.  (30),  but  with  x being 
replaced  by  the  arc  length  coordinate,  s . For  the  case  of  a laminar-flow 
body  we  assume  a sine  function  velocity  profile,  which  according  to  Schlichting 
[6],  very  accurately  approximates  the  exact  Blasius  profile.  Thus,  we  assume 


J1  " “V 


tt  (r  - rfe) 


6.  = 4.7953  (vs/U  ) ' . 
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With  these  assumed  velocity  profiles,  Eq.  (32)  can  be  integrated.  Without 
showing  the  details,  we  find: 
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C1  “ uv  1 
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1 2 16  9 w w 


2 2 
r - r. 
w b 
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for  laminar  flow  on  the  body,  and 


[2  " "v 


for  turbulent  flow  on  the  body 
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We  now  eliminate  Uy  from  1^  ; thus,  Eq.  (12)  must  be  integrated.  Again, 
we  write 


V*i 
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w w w 


*i  = | u rdr  =|  u^  rdr  + Uy  j rdr  + uw  rdr  , 
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Because  * Aq  U^^tt  we  obtain  for  i * 1 (laminar  flow): 

"v  " *«  "■  I » - I I 5 2 + 2(2.„)  r 5 . + ■ { 2 

1 b 1 4 15  w 


* 

f \ 
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and  for  i = 2 (turbulent  flow) 
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The  first  term  of  ACp  may  be  expressed  in  the  form  ACp^=ACp^i(i=l  or  2). 
For  laminar  flow  on  the  body,  we  find: 
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and  for  turbulent  flow  on  the  body,  we  find: 
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If  the  boundary  layers  are  thin,  one  may  of  course  neglect  the  higher  order 


terms  in  6 and  <5,  . 
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APPLICATION  OF  THE  THEORY 

Body  Shape 

Based  on  the  above  analysis,  computations  have  been  performed  for  a large, 
streamlined  body  designed  to  operate  in  the  Garfield  Thomas  1,22— m.  diameter 
Water  Tunnel  [8].  This  body  was  referred  to  in  the  INTRODUCTION  and  is 
described  mathematically  by: 


y » /x(2-x)  - K 


C x 
o 

n 2a2 


exp< 


~2 
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2a 


2-kKn 


+ x2  eY 


(46) 


where  x = x/A  , A - nose  length,  y = r./D_/2,  D - maximum  diameter  of  the 
n n b max  max 

bodv.  K = A /D  . and  a,  C , and  k are  numerical  constants.  The  parameter 
n n max  o 


e 


L 


K C 
n o 

2a  ^ 


expi 


(47) 


The  tailcone  is  formed  by  an  18°  arc  of  radius 


2.5  D 


*L 


max 


2 sin  18° 


(48) 


An  inflection  curve  is  then  faired  between  the  end  of  the  arc  (at  4>=90— 18=7 2 ) 

and  a 4-inch  diameter  disk  whose  center  is  coincident  with  the  axis  of  symmetry. 

This  body  shape  was  shown  in  Figure  1 for  a = 0.3,  Cq  = 0.0303,  and,  k = 0.45227. 

The  total  length,  L , is  3.05-m.  and  D is  32.39  cm. 

* max 


Experimental  Pressure  Distributions 

The  experimental  body  was  manufactured  by  building  up  multiple  layers  of 
fiberglass  cloth  and  resin.  A precision  template  was  used  to  screed  the  body 
to  final  shape  before  the  last  coating  of  resin  fully  hardened.  Figure  3 shows 
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a photograph  of  the  completed  body  and  sting  mount  used  to  hold  it  in  the  1.22-m. 
diameter  test  section.  A total  of  94  pressure  taps  were  located  in  the  surface 
of  the  body.  These  pressure  orifices  were  0.  79-mm.  in  diameter.  Each  body 
tap  was  connected  to  a flexible  plastic  tube  which  passed  through  the  inside  of 
the  body,  through  the  sting,  and  out  of  the  tunnel.  The  94  tubes  then  connected 
to  two  48-channel  Scanivalves  and  a single  Bell  and  Howell  calibrated  pressure 
transducer.  Before  all  runs,  the  tubes  were  bled  to  remove  air  pockets.  Figure  4 
shows  a schematic  of  the  experimental  arrangement  and  the  equations  used  in  data 
reduction  to  obtain  Cp  . 

£i 

As  indicated  in  Figures  2(c)  and  4 the  reference  velocity  in  the  measurements, 
u^,  was  obtained  with  a pitot  probe  placed  relatively  close  to  the  body  and  tunnel 
wall.  Consequently,  u^  may  not  be  interpreted  as  the  used  in  the  theoret- 

ical analysis.  The  velocity  u^  differs  from  Uw  because  of  both  potential 
effects  induced  by  the  body  and  viscous  effects  created  by  tunnel  wall  boundary- 
layer  growth. 

Detailed  measurements  of  the  turbulent  boundary  layer  growth  on  the  wall  of 
the  Garfield  Thomas  Water  Tunnel  were  performed  and  reported  by  Ross,  et  al  [7]. 
They  concluded  that  an  effective  origin  for  the  start  of  boundary  layer  growth  is 
1.5  rw  upstream  from  the  geometric  beginning  of  the  test  section.  That  is,  one 
may  assume  that  the  cylindrical  test  section  can  be  ficticiously  extended  up- 
stream into  the  settling  section  a distance  1.5  rw  and  refer  to  this  point  as 
the  entrance  of  the  test  section.  Figure  2(c).  Thus,  we  can  assume  £ **1.5r  +d  , 
where  dQ  is  the  distance  from  the  nose  of  the  body  to  the  actual  beginning  of 
the  cylindrical  test  section  (dQ  * 54.6  cm.  in  this  arrangement). 

In  order  to  compare  our  experimental  pressure  distributions  with  theoreti- 
cally predicted  ones,  it  is  apparent  that  Cp_  (referenced  to  u ) must  be 

h m 


GCL:mmJ 


corrected  in  order  to  account  for  the  differences  between  u and  U . We 

m 00 

do  this  in  two  steps. 

We  first  correct  Cp  for  the  tunnel  wall  boundary-layer  displacement 
thickness  at  the  position  of  the  pitot  probe.  This  thickness  varies  from 
2.02-mm.  at  15.24  m./sec.  to  2.42-mm.  at  6. 09m. /sec.  (the  range  of 
velocities  considered  in  this  investigation).  The  pitot  probe  was  4-inches 
from  the  wall  and  thus  was  out  of  the  boundary  layer.  The  effective  area  of 
the  tunnel  at  the  position  of  the  pitot  probe  is: 


A - 2tt  r 6 + tt  6 ^ 

o w w w 


= A - 2tt  r 6 
o w w 


(49) 


From  the  continuity  equation  [Eq.  (12)],  we  can  write 
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(50) 


which  represents  the  increase  in  relative  to  due  to  boundary-layer 

growth.  The  experimental  pressure  coefficient  corrected  to  Uw  by  viscous 
considerations  only  is  thus: 
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We  know  (because  of  the  close  proximity  of  the  pitot  probe  to  the  body) 
that  u^  may  also  be  influenced  by  the  potential  effects  of  the  body.  In 
order  to  correct  for  this,  a potential  flow  solution,  such  as  that  obtained 
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from  the  Douglas-Neumann  code  is  required.  We  calculate  the  velocity  ratio 

u = U /u  at  the  pitot  probe  location.  Then, 
rp  00  m ’ 


Cp'°r-  - 1 - a - Cp“)/U^  (52) 

is  the  final  correction  of  Cp^  . Furthermore  it  makes  no  difference  which 
correction  is  applied  first,  i.e., 

iet  £Pe  - 1 - <1  - CPE)/ur2p 

then  Cpg°r‘  = 1 - (1  - Cp*) 
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e rp 
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It  is  easily  verified  that  Eq.  (52)  also  reduces  to  Eq.  (53).  For  the  experi- 
mental arrangement  used  here,  we  found  u = 0.9976,  and  that  A /A  ranges  from 

rp  o e 

1.013  to  1.016. 

The  experimental  pressure  distributions  shown  in  Figure  1 have  been 
corrected  according  to  the  above  discussion. 


Results 

The  generality  of  the  above  analysis  allows  both  laminar  and  turbulent- 
body  flow.  As  part  of  the  experimental  program,  we  therefore  measured  the 
laminar-to-turbulent  transition  points  at  the  four  tunnel  velocities  in  which 
pressure  distributions  were  measured.  A wedge-shaped  hot-film  probe  and 
anemometer  system  were  used  in  these  measurements.  TABLE  I summarizes  the 
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transition  location  results.  For  those  body  coordinates  forward  of  st 
(transition  arc  length),  the  laminar  form  of  the  pressure  coefficient 
correction  analysis  is  used  [Eq.  (44)].  And,  in  Eq.  (26)  is  calcu- 

lated according  to  the  formula 

Tb  * 0.332  pvUa  (U J\)  s)1/2  (54) 

in  this  laminar-flow  regime.  For  s > st»  is  calculated  according  to 

Eq.  (29)  with  the  appropriate  change  from  x to  s variable.  The  integral 

of  Eq.  (26)  is  evaluated  using  Simpson's  Rule  of  numerical  integration. 

The  integrals  of  Eq.  (28)  are  also  evaluated  using  Simpson's  Rule. 

Both  pb  and  p^  are  determined  from  the  body-in-a  cylindrical  duct  Douglas- 

Neumann  solution.  In  this  numerical  example  we  calculated  <5*  (needed  to 

b 

obtain  p^)  using  an  axisymmetric,  implicit  finite-difference  solution  of  the 
boundary  layer  equations  of  motion. 

Because  of  the  thicker  boundary  layers,  ACp  is  largest  at  the  lower 


velocities.  Thus,  at  6.09  m./sec.  a few  typical  values  of  ACp^(i=l,2,3,4) 
are  presented  in  TABLE  II  from  which  an  assessment  of  each  ACp^'s  contribu- 
tion to  the  total  ACp  can  be  made.  Clearly,  ACp^  contributes  the  greatest 
to  ACp  , but  ACp^  and  ACP2  are  also  significant  contributors.  On  the 
other  hand,  ACp^  is  approximately  an  order  of  magnitude  lower  than 
ACp2(ACp2  < ACp^  < ACp^)  for  those  locations  on  the  body  where  the  flow  is 
turbulent.  Obviously,  ACp^  is  the  most  tedious  to  calculate,  so  it  may  be 
conveniently  justified  to  neglect  it  entirely  from  the  computations.  For 
those  conditions  in  which  the  computations  are  compared  to  experimental  data, 
this  deletion  has  been  made.  It  is  further  noted,  that  if  body-flow  is 
entirely  laminar,  ACP2  may  also  be  neglected;  although,  we  have  retained 
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it  for  all  of  the  computations  in  this  example  because  there  is  only  a small 
regime  of  laminar  flow  on  the  body. 

Figures  5 through  8 summarize  the  comparisons  of  potential  flow,  corrected 
potential  flow,  and  experimental  pressure  distributions  for  our  given  body 
operating  at  6.09  through  15.24  m. /sec.  test-section  velocities.  At  all  speeds 
the  horizontal  buoyancy  correction  theory  agrees  well  with  the  experimental 
data,  particularly  for  those  body  positions  forward  of  CpM<n  . Aft  of  Cpm^n , 
there  is  a region  of  strong  interaction  where  coupling  between  the  boundary- 
layer  flow  and  the  potential  outer  flow  becomes  important.  The  theory  presented 
here  is  not  applicable  to  this  region. 

CONCLUSIONS 

The  analysis  presented  here  is  applicable  to  many  practical  problems 
involving  the  flow  about  bodies  of  revolution  within  cylindrical  ducts,  such 
as  wind  or  water  tunnel  test  sections.  In  particular,  this  analysis  allows 
one  to  correct  an  inviscid,  in-duct  body  pressure  distribution  for  the  effects 
of  horizontal  buoyancy.  The  derived  expressions  are  relatively  simple;  thus, 
computations  are  performed  using  very  simple  Fortran  or  desk-top  calculator 
programs.  Based  on  the  experimental  observations  of,  and  on  the  computations 
performed  for  a particular  body  that  operates  in  a given  water  tunnel  test 
section,  horizontal  buoyancy  effects  were  shown  to  be  important  and  should 
not  be  neglected  if  precise  pressure  coefficient  values  are  required.  We  have 
demonstrated  that  the  theory  adequately  predicts  the  observed  effects,  except 
in  the  afterbody  region  where  strong  interaction  effects  dominate. 
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Table  I The  experimentally  determined  transition  points  on  the  body  used  in 
the  pressure  distribution  measurements  . 


Tunnel  Velocity 
U«, 

(m./sec.) 

Body-Length  Reynolds  No. 

Re.  = U L/V 

L 00 

6.09 

1.88  x 107 

9.14 

2.83  x 107 

12.19 

3.76  x 107 

15.24 

4.72  x 107 

Transition  Point 


0.083 


Transition 
Reynolds  No 
ReT (s_/L) 


2.82  x 10 
3.28  x 10 
3.76  x 10 
3.93  x 10 


Table  II  Typical  values  for  ACpj  (i  = 1,  2,  3,  4)  at  6.09  m./sec., 
inviscid  pressure  coefficient,  and  the  corrected  inviscid 
pressure  coefficient  for  horizontal  buoyance  effects. 


Body  Coordinate 

x/L  = (x-Jl  )/L 
o 


ACp2  ACp3  ACp4  ACp  = 1 ACp±  Cpj 

i=l 


.002  -0.0038  0.000  -0.0139  0.00009 

.075  -0.0046  -0.00003  -0.0158  0.00047 

.253  -0.0065  -0.0006  -0.0202  0.00061 

.497  -0.0091  -0.0020  -0.0261  0.00069 

0.667  -0.0108  -0.0031  -0.0300  0.00071 


-0.0177 
-0.0199 
-0.0267  1-0.087 


-0.0366 

-0.0432 


-0.139 

-0.161 


PRESSURE  COEFFICIENT 
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Figure  1 The  pressure  distributions  on  a large  axisynunetric  body 

(L  * 3.05-m.)  operating  in  a cylindrical  test  section 

(r  » 1.22-m.) 
w 


NONDIMENS IONAL  RADIUS  (r/r^ 


Figure  2 (a)  Potential  flow  over  a body  in  an  infinitely  long  cylindrical  duct, 

(b)  a schematic  description  of  boundary  layer  buildup  on  the  body  and 
duct  wall  depicting  that  uv  is  different  from  u^  and  (c)  a schematic 
description  of  the  flow  over  a body  in  an  experimental  situation. 
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Figure  4 A schematic  of  the  data  acquisition  and  reduction  requirements 
used  to  obtain  the  experimental  pressure  distributions. 
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Figure  5 A comparison  of  the  in-tunnel  potential  flow  and  corrected 
potential  flow  pressure  distributions  with  one  determined 
experimentally  at  6.09  m./sec.  in  the  1.22-m.  diameter 
water  tunnel. 


NONDIMENS fONAL  RADIUS  <r IrJ 


PRESSURE  COEFFICIENT 


-33- 


June  1,  1978 

GCLrmmj 


Figure  6 A comparison  of  the  in-tunnel  potential  flow  and  corrected 
potential  flow  pressure  distributions  with  one  determined 
experimentally  at  9.14  m./sec.  in  the  1.22-m.  diameter 
water  tunnel . 
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Figure  7 A comparison  of  the  in-tunnel  potential  flow  and  corrected 
potential  flow  pressure  distributions  with  one  determined 
experimentally  at  12.19  m./sec.  in  the  1.22-m.  diameter 
water  tunnel. 
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Figure  8 A comparison  of  the  in-tunnel  potential  flow  and  corrected 
potential  flow  pressure  distributions  with  one  determined 
experimentally  at  15.24  m./sec.  in  the  1.22-m.  diameter 
water  tunnel. 
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